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Verfahren zur Bewertung der Betriebsbedingungen einer 
Maschine oder einer Anlage 

Die Erfindung betrifft ein Verfahren zur Bewertung der 
5 Betriebsbedingungen einer Maschine oder einer Anlage . 

Bisher wurden zur Bewertung der Betriebsbedingungen einer 
Maschine oder einer Anlage visuelle Extrapolationen eines, 
bspw. kritischen, Messwertverlauf s und/oder subjektive 

10 Bewertung en der Auswirkungen auf Parameter der Maschine oder 
der Anlage vorgenommen, urn Betriebsbedingungen einer Maschine 
oder einer Anlage abzuschatzen und entsprechend zu reagieren, 
indem der Parameter entsprechend der Bewertung bspw. durch 
Vorgabe eines Sollwertes fur Parameter verandert oder in 

15 sons tiger Art und Weise beeinflusst wird. 

Dabei wurden Anpassungs-/Regressionsfunktionen einer 
Datenbasis angepasst und eine Optimierung durch iterative 
Auswahl von Kurvenf unktionen, welche den maximalen 
Korrelationskoef f izienten auf weist , durchgef iihrt . 

20 Der Kriimmungsverlauf einer derart angepassten Kurve muss 
nicht zwingend mit dem der Datenbasis tibereinstimmen. 
Der Korrelationskoef fizient r (maximaler Wert = 1, minimaler 
Wert = - 1) ist als Anpassungsgutekriterium nur bedingt 
brauchbar, da dieser Wert nicht nur von der Anpassungsgute 

25 der angewendeten Kurvenf unktion abhangt, sondern auch von der 
Steigung der verwendeten Kurvenfunktion . Lauft die Steigung 
z.B. einer linearen Anpassung gegen null, folgt dieser 
Tendenz auch r, unabhangig von der Streuung der einzelnen 
Kurvenpunkte . 

30 Somit ist r als GiitemaS fur eine Extrapolation nicht 
brauchbar . 

Werkzeuge zur Generierung fur die Bewertung der 
Betriebsbedingungen einer Maschine oder Anlage mlissen 
35 bestimmte Mindestanf orderungen erfullen, damit eine Bewertung 
des Extrapolationsergebnisses hinsichtlich der 
• Vorhersagezuverlassigkeit 



200206132 



• Einf lussgroSen auf eine Vorhersagesicherheit 

• Nachvollziehbarkeit der Vorhersage 
moglich ist. 

5 Naturgemafi ist jede Vorhersage mit Unsicherheiten behaftet 
und deren GutemaS zur Beurteilung von prognostizierten 
Alarmen/Grenzwertiiberschreitungen sowie als 
Entscheidungsbasis fur daraus resultierende bspw. 
automatisierte Aktionen auSerordentlich wichtig. 

10 Dariiber hinaus kann bei einem zyklisch generierten Verfahren 
zur Bewertung der Betriebsbedingungen einer Maschine oder 
Anlage eines bestimmten, z.B. messtechnischen Zusainmenhangs 
die zeitliche Entwicklung eines GiitemaSes wiederum als Trend 
betrachtet werden und damit zusatzliche Aufschlusse uber den 

15 Eintref f zeitpunkt des prognostizierten Ereignisses, z.B. 
einer Grenzwertiiberschreitung, geben, urn eine 
Betriebssicherheit der Maschine oder der Anlage zu 
gewahrleisten . 

Die Einschatzung der Zuverlassigkeit einer Extrapolation 
20 besitzt zentrale Bedeutung, da Bewertungen von 

Schadenssymptomen, Grenzwertiiberschreitungen und 
Betriebsoptimierungen, z.B. Ersatz von Verbrauchsteilen, 
ohne Kenntnis der Vorhersagezuverlassigkeit , also des 
GutemaSes kaum moglich sind. 

25 

Eine Bewertung der Betriebsbedingungen einer Maschine oder 
einer Anlage, bei dem ein Trend der weiteren Entwicklung von 
Parametern analysiert wird, ist sehr wichtig im Zusammenhang 
mit Problemlosungen bei Auf gabenstellungen wie 

30 Betriebsmonitoring oder Systemanalyse . 

Diese Verfahren werden uberall dort eingesetzt, wo 
Extrapolationen in unbekannte Wertebereiche aus bekannten 
GesetzmaSigkeiten heraus von Bedeutung sind, also im 
wesentlichen bei Ereignisprediktionen (Erkennung zukiinf tiger 

3 5 Alarme/Uberschreitung von Belastungsgrenzen) und 

Schadenspraventionen (Erkennung zukiinf tiger Schaden) . 
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Ahnlich wie bei Messwerten, diese sind ohne Angabe der 
Messtoleranz praktisch wertlos, muss einer Extrapolation 
eine GroSe zugeordnet werden, welche Aufschluss iiber die 
5 Vorhersagequalitat (GlitemaS) gibt. Ohne diese Information 
sind Extrapolationsergebnisse kaum bewertbar und daher 
Entscheidungen zu Folgeaktionen (z.B. Prozesseingrif f e) 
hochst unsicher und u.U. kontraproduktiv. 

10 

Es ist daher Aufgabe der Erfindung dieses Problem zu beheben. 



Die Aufgabe wird durch ein Verfahren gemass Anspruch 1 
15 geldst. 

Weitere vorteilhafte Ausgestaltungen sind in den 
Unteranspriichen aufgelistet. 

20 

Die optimale Anpassung der Extrapolationsf unktion an die 
Datenbasis erfolgt nicht durch einen maximalen 
Korrelationskoef f izienten r, sondern mit Hilfe eines 
GutemaSes K, welches von der Steigung einer Anpassungskurve 
25 an eine Datenbasis unabhangig ist. 

Aus mindestens zwei verschiedenen Unterkriterien wird das 
GutemaS K berechnet, welches Aufschluss iiber die 
Vorhersagequalitat der Extrapolation erteilt. 
3 0 Die Existenz von K ermoglicht zumindest zwei 
Extrapolationsmodi : 

• Vorgabe eines konstanten GiitemaSes K 

3 5 In diesem Falle ist die Extrapolationsweite, d.h. der 

Abstand zwischen dem letzten Wert der Datenbasis und einem 
Endpunkt der Extrapolation, variabel . 



4 

Diese Vorgehensweise ist wichtig bei der Unterdrttckung von 
Fehl- oder unbedeutenden Alarmen. 

Hier kann einem Operateur einer technischen Anlage 
mitgeteilt werden, iiber welchen Zeitraum kein Alarm 
auftreten wird, d.h. eine bestimmte Betriebssicherheit 
gegeben ist. 

Vorgabe einer Extrapolationsweite (z.B. bei 
Grenzwertuberschreitung ) 

in diesem Falle ist K variabel . Die zeitliche Entwicklung 
von K kann wiederum als Basis fur eine Trendanalyse 
benutzt werden, um den wahrscheinlichen Eintref f zeitpunkt 
des Ereignisses bei K = 1 abzuschatzen. 



Trendanalysen, um die Betriebsbedingungen einer Maschine oder 
Anlage zu bewerten, konnen ganz allgemein die Entwicklung 
zweier korrelierter GroSen in noch unbekannte Werteberexche 
extrapolieren . 

Ihre Anwendung kann sowohl auf zeitliche Entwicklungen 
yi = f (ti) ; i = 1. . .n 

Beispiel: P e i=f(t); P el = Generatorlei stung 

und auf formal zeitunabhangige Zuordnungen 

yii = f(y 2 i); 1 = 1 --- n ' yi mit Y2 synchron 

Beispiel: P e l = f (Brennstof fmassenstrom) 
angewendet werden. 

Der Einsatz von Trendanalysen unterstellt allerdings 
implizit, dass sich die aus der Vergangenheit ableitbare 
GesetzmaSigkeit der Entwicklung zweier korrelierter GroSen 
auch in der Zukunft fortsetzt und sich auf diese Weise noch 
unbekannte Wertezuordnungen generieren lassen. 
Diese Bedingung setzt far sinnvolle Trendanalysen voraus, 
dass diese nur mit Verlaufen generiert werden konnen, welche 
unter identischen ( ,normalisiertenM Randbedingungen - 
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physikalische Umweltbedingungen, Brennstoff, Leistung, 
Betriebsmodus Brenner usw. - generiert wurden. 
Die Bereitstellung solch normal isierter Wertezuordnungen kann 
entweder durch Datenselektion mit Hilfe logischer Bedingungen 
5 (Diskrete Normalisierung) oder mit Hilfe eines 

Prozesssimulationsmodells (analytische Normalisierung) 
erfolgen, welches z. B. alle umgebungsabhangigen Werte auf 
ISO-Bedingungen zuriickrechnen kann. 

10 

Es zeigen: 

Figur 1 eine Datenbasis und Extrapolationskurve, 
Figur 2 Kurven zur Bestimmung von Al 
Figur 3 eine Segmentierung einer Datenbasis, 
15 Figur 4 einen zeitlichen Verlauf von K, 
Figur 5 eine Turbine mit Bauteilen, 

Figur 6 in vereinf achter Darstellung einen Langsschnitt einer 
Turbine , und 

Figur 7 eine Turbinenschauf el als ein Bauteil. 

20 

BasisgroSen zur Berechnung zum Gutemafi K 

Das GutemaS K ist ein Ma£ fur die Giite einer Extrapolation, 
2 5 Sie verkmipft mehrere Fehlereinf lusse und ist daher eine 
Funktion zumindest zweier Variabler z. B. V, Al . 
In diesem Falle sind dies bspw. ein Verhaltnis V, eine x- 
Unscharfe Al, eine Stetigkeit S und eine zeitliche Konstanz 
C: 

30 

K = f (V, Al, S, C) . 

Diese GroSen werden bspw. so verkniipft, dass sich der 
Wertevorrat von K innerhalb eines normierten Bereiches, z.B. 
35 0% bis 100% befindet. 
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• Verhaltnis V des Datenbereiches zum Extrapolationsbereich 

Fur die Extrapolation stehen n Messpunkte (x^yi) bis (x n ,y n ) 
zur Verfiigung. Es miissen jedoch nicht alle vorhandenen 
5 Messpunkte 15 eines Datenvorrats 54 (Fig. 3) verwendet 
werden . 

Das Verhaltnis V setzt bspw. den x-Wertebereich (x n -xi) der 

zur Extrapolation verwendeten Datenbasis 6 in das Verhaltnis 

der Strecke xi bis zur x-Position des verwendeten 
10 Extrapolationspunktes x s (Fig.l), d.h. dem 

Extrapolationsbereich 9, in ein Verhaltnis, welches bspw. 
maximal gegen den Wert 1, d.h. die Datenbasis 6 ist gleich 
dem Extrapolationsbereich, und minimal gegen den Wert 0, d.h. 
die Datenbasis 6 ist sehr klein gegeniiber dem 
15 Extrapolationsintervall , strebt . 

X s ist z.B. ein Wert, bei der eine Grenzwertiiberschreitung 
(Y s = Grenzwert 18) stattfindet. 

In Figur 1 ist ebenfalls eine Anpassungsf unktion 12 an die 

Datenbasis 6 gezeigt, die auch die Extrapolation bis zum x- 
20 Punkt x s ergibt. 

V wird bspw. wie folgt berechnet: 

V = (x n - xi)/(x s - xi) . 

25 

Tragt V wesentlich zu einer Reduzierung des GutemaSes K bei, 
kann schrittweise die Datenbasis der Extrapolation erweitert 
werden . 

30 • Unscharfe Al der Anpassungskurve in x-Richtung 

Jede Kurvenanpassung (Voraussetzung zur Extrapolation) 
besitzt naturgemafi auch eine Unscharfe (Vertrauensbereich) 
der berechneten Kurvenparameter . 
35 Urn die Unscharfe der Steigung zu bestimmen, wird in der Figur 
2 bspw. eine lineare Regressions f unktion 21 verwendet. 
Generell konnen als Extrapolations funktionen alle Funktionen 
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verwendet werden, welche sich in lineare Strukturen 
uberfuhren lassen, also 

Lineare Funktion (Stammfunktion) -> y = ap + a^*x 

Potenzfunktion -> In y = lnag + a]_*lnx 

5 Logarithmische Funktion -> y = ag + a]_*ln x 

Exponentielle Funktion -> In y = In ag + a x *x 

u.a. . 

Die UnschcLrfe Al wird dann wie folgt berechnet. 
10 Figur 2 zeigt eine beispielhaf te Datenbasis 6 mit einer 

Extrapolationskurve 21 und mit weiteren Kurven 24, 27 zur 
Berechnung von Al. 

Die Unscharfe Al kann durch Drehung der Anpassungskurve 21 
15 (in linearer Darstellung) urn einen Drehpunkt (x,y) 

dargestellt werden. Die Mittelwerte jcund y sind bspw. durch 
geometrische oder arithmetische Mittel bestimmt worden. Der 
Punkt (x y y) muss nicht zwingend auf der linearen 

Anpassungskurve 21 (Extrapolationskurve) liegen. 
20 Die Unscharfe einer Steigung der linearen Kurve 21 entspricht 
einer Winkeldrehung, bei der die Steigung m urn ±Am variiert 
wird und damit einer x-Unscharfe Al der Anpassungskurve u.a. 
in x-Richtung (Fig. 2) . 

Durch die Variation der Steigungen ergeben sich zwei weitere 
25 Kurven 24, 27, die mit einer Parallelen 18 (ys=konstant) , die 
einem Grenzwert entspricht, jeweils einen Schnittpunkt 30, 33 
bilden. Zu jedem Schnittpunkt 30, 33 gibt es einen 
entsprechenden x-Wert Imax und I m in/ wobei Ima X > I m i n ist und 

Al — Imax ~~ Imin • 

30 Diese x-Unscharfe korreliert mit einem Schnittwinkel a, der 
sich bestimmt aus dem Winkel zwischen der Parallelen 18 und 
der Kurve am Schnittpunkt der Extrapolationskurve 21 und der 
Parallelen 18. Die x-Unscharfe Al verstarkt sich mit 
abnehmendem Schnittwinkel a. 

3 5 Ist bei einer Extrapolation kein Grenzwert vorhanden bzw. 
wird keiner uberschritten, gilt das oben gesagte fur den 
vorgegebenen x-Endpunkt (x s ) der Extrapolation und einer 
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gedachten Parallelen zur x-Achse, welche durch den Punkt 
(x s ,y s ) lauft. 

Da die Extrapolat ions funkt ion immer analytisch vorliegt, kann 
mit Hilfe der ersten Ableitung der Schnittwinkel a zwischen 
5 einer Waagerechten oder Alarmfunktion und der 
Extrapolationsfunktion berechnet werden. 

Der Ansatz zur Berechnung der Steigungsunscharf e ±Am ist 
bspw. wie folgt: 
10 Eine Steigung b wird mittels eines Konf idenzintervalls 
bestimmt. Ein solches Verfahren ist aus Kreyszig, Erwin: 
"Statistische Methoden und ihre Anwendungen" , Verlag 
Vandenhoeck und Ruprecht, Gottingen, Seite 270 bekannt. 

15 Zur Bestimmung der Steigung b wird wie folgt vorgegangen: 

a) Eine Datenbasis 6 wird festgelegt. Die Datenbasis 
6 umfasst n korrelierte x- und y-Werte (Fig. 2) 

20 b) Berechnung von jc und y der Datenbasis 6 sowie 

der Grofie ^x.y t 

c) Berechnung von =— — ( ^x^ -nxy) (i=l..n) 

d) Berechnung von e i- — — ( Yx - x) 2 (i=l..n) 

x n-l 

e) Berechnung von ez = — ( Yy l -y) a (i=l..n) 

y n-l 

25 f)die Steigung b ergibt sich aus b = -^y 

S x 

g)Berechnung von a = (n-l ) ( -b 2 s 2 x ) 

Somit ergibt sich die Gleichung der Regressionsgeraden 21 
30 y = y + b(x-x) . 

Als nachstes wird eine Konf idenzzahl von bspw. 95% 
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vorgegeben, d.h. y = 0,95, woraus eine Grosse F (c) berechnet 
wird: 

h) F(c) = ^(l + r) = 0,975 

i) mit F(c) = 0,975 und n-2 (n = Anzahl Messwerte) 
5 Freiheitsgrade liefert die t-Verteilung (Student- 

Verteilung) einen Wert c 

(0,975 entspricht dem Integral der t-Verteilung bis 
an die Stelle x=c) . 

10 Am ergibt sich aus 
. j c4a 

Somit ergibt sich eine Unscharfe der Steigung m: 

15 b - Am < m < b + Am . 

Berechnet wird damit die x-Unscharfe Al der Steigung m im 
Drehpunkt (Fig. 2) in Abhangigkeit der StreuungskenngroSen 
der aktuellen Datenbasis 6. 
20 Damit ist Al = f (Am, a) mittels Strahlensatz oder 

Goniometrie berechenbar . 

Eine Normierung von Al kann ebenso erfolgen. 

25 • Stetigkeit S der y-Werte der Datenbasis 

Eine Segmentierung einer Datenbasis 6 erteilt vorab 
Aufschluss iiber die aktuelle Kriimmungsrichtung der 
verwendeten Datenbasis 6 . 
30 Die Datenbasis 6 ist eine echte oder unechte Teilmenge 

eines Datenvorrats 54, der alle verfiigbaren Messpunkte 15 
umf asst . 

Eine Segmentierung bedeutet, dass Messpunkte, also 
Zahlentupel (x,y) der Datenbasis 6 in zumindest drei 
3 5 Segmente 45, 48, 52 aufgeteilt werden. Fur jedes Segment 
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44, 48, 52 wird eine lineare Anpassungskurve 36, 39, 42 
bestiirant . 

Die Kriimmungsrichtung der Datenbasis 6 wird durch Bildung 
von Segmentmittelwerten gl, g2 und g3 und Berechnung der 
5 zweiten numerischen Ableitung bestiirant. 

Bei diesem Verfahren kann davon ausgegangen werden, dass 
die Berechnung der Anpassungskurve (Extrapolationskurve) 
zyklisch aus einer gleitenden Datenbasis 6 erfolgt, d.h. 
10 ein Datenfenster 3 (Fig. 3) konstanter oder variabler 

Lange wird nach jedem abgeschlossenen Extrapolations zyklus 
in Richtung zunehmender x-Grofien um ein jeweils 
vorgegebenes Interval 1 Ax (unabhangige Variable) 
verschoben. 

15 In der Praxis kann der Vorschub Ax des Datenf ensters 3, 

d.h. der Datenbasis 6, maximal in dem zeitlichen Zyklus 
erfolgen, mit welchem neue Messwerte 15 generiert werden 
( z . B . Ax = 5 s ) . 

Der Begriff der Stetigkeit S ist in diesem Zusammenhang 
20 nicht als mathematische Definition zu verstehen, sondern 

als das Mafi der Veranderung des Punktverlauf es der 
Datenbasis 6 in Relation zum letzten abgeschlossenen 
Extrapolationsschritt . 

Um eine Aussage iiber die Stetigkeit der Y-Werteentwicklung 

25 und damit implizit eine Aussage iiber die 

Extrapolationsfahigkeit der Entwicklung machen zu konnen, 
wird das Datenfenster 3 (Fig. 3) in mindestens drei 
Segmente 45, 48, 52 unterteilt. Fur jedes Segment 45, 48, 
52 wird in einem Zwischenschritt je ein Mittelwert und je 

30 eine lineare Anpassungsfunktion yi, y 2 , y3 mit den 

Steigungen c 1# c 2 und c 3 berechnet (Fig. 3) . 
Bezeichnet man die Mittelwerte der drei Regressionsgeraden 
36, 39, 42 mit g lf g 2 (mittleres Segment) und g 3 , so kann 
mit dem numerischen KrummungsmaS p 

35 P = gi - 2*g 2 + g 3 

die aktuelle Kriimmungsrichtung des Steigungsverlauf es 
bestimmt werden. 
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Im Beispiel von Figur 3 ist q negativ. Dies bedeutet eine 
Rechtskrummung . 

Bei p = 0 liegen die drei Mittelwerte auf einer Geraden. 
Aus einem Kurvenvorrat wird iterativ jener Kurventyp dem 
5 Wertebereich der gesamten aktuellen Datenbasis angepasst. 

Anpassungs funktionen (Regress ions funktionen) des 

Kurvenvorrats fur die Extrapolationskurve von 

Datenkorrelationen miissen die Bedingung des monotonen 

10 Verlaufes erftillen, da nicht -monotone Funktionen zu 

erheblichen Extrapolationsunsicherheiten fuhren konnen. 

Generell sind alle Funktionen verwendbar, welche in 
lineare Strukturen der Form y = ag + a^*x uberfiihrbar 

sind, also 

15 

Lineare Funktion(Stammfunktion) -> f(x)=y = ag + a]_*x 

Potenzfunktion -> f(x)=lny= lnag + ai*lnx 

Logarithmische Funktion -> f (x)=y = ag + a^*lnx 
Exponentielle Funktion f(x)=lny = In ag + al*x 

20 u.a. , 

Der aus dem Kurvenvorrat ausgewahlte Kurventyp muss 
folgende Bedingungen erfiillen: 

25 die Krummungsrichtung der Extrapolationskurve 

(Regressionsgeraden) muss mit der von p ubereinstimmen 

der Quotient aus Zahler = (ggf . gewichteter) Mittelwert 
der Abstandsguadrate zwischen Messwerten und 
30 Extrapolationskurve, und Nenner = Mittelwertquadrat des y- 

Wertebereiches der Extrapolationskurve im Bereich des 
Datenf ensters (dabei ist beispielsweise y 2 mitt k = t (y maXK 
+ y m in K )/2] 2 , wobei ymax der maximale Y-Wert ist und ymin 

entsprechend der minimale Y-Wert der k-ten Kurve) und 
35 dient zur Normierung) muss minimal sein: 
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Q k = f(£) = — : r^e = min (i = l.-nim) 

y 2 . * > w. 

^ mm, » 



k ist dabei ein Zahler der verfiigbaren 

Extrapolationskurventypen (Kurvenvorrat ) . 
5 Yi(xi) ist der Messwert an der Stelle x-^ . 

fk(x^) ist der Funktionswert des k-ten 

Extrapolationskurventyps an der Stelle x^. 

W£ ist ein Wichtungsf aktor fiir jeden einzelnen Messwert oder 
fur alle Messwerte eines Segments. 

10 

Die Bewertung der Stetigkeit S wird wie folgt vorgenommen: 
Vergleich der drei Geradensteigungen Ci, C 2 , C3 mit den 
Steigungen der Extrapolationskurve an den jeweiligen 
Mittenpositionen (X s i,X s2 , X s3 ) der drei Segmente des 
15 Datenf ensters . Eine unterschiedliche Wichtung der drei 

Steigungsdif ferenzen ist hier moglich und sinnvoll. So konnen 
die aktuelleren Werte im letzten Segment starker bewertet 
werden, um schneller eine Veranderung des Kurvenverlauf s zu 
erkennen . 

20 Die drei Steigungsdif ferenzen sind ein MaE fiir die 
Entwicklungsstetigkeit des Kurvenverlauf es . 



S berechnet sich wie folgt: 

Oi bis O3 sind die drei Steigungen der ausgewahlten k-ten 
25 Anpassungskurve 36, 39, 42 an den jeweils halben 
Segment sbr e i t en , 

Ci bis C3 die Steigungen der linearen Segmentsanpassungen 
(Fig. 3) . 



30 Somit gilt: 



Yyi*{C x -o x y 

S = — — =; ; i = 1...3 (Segmentbereiche) 

2>< 



mit Yi : Wichtungsf aktoren l...n. 
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Der Wertebereich von S ist wie folgt: 0 (abs. Stetigkeit) , 
d.h. die Steigungen der Extrapolations funkt ion in den 
Segment smitten sind identisch mit den Regressionsgeraden der 
5 einzelnen Segmente 45, 48, 52 bis + «>, d.h. keine 
Ubereinstimmung der Steigungen. 

• Zeitliche Konstanz C der Extrapolationen 

10 

Die zeitliche Varianz einer Folge von Extrapolationen ist ein 

zusatzlicher und wichtiger Indikator der 

Vorhersagesicherheit , da hohe Varianzbetrage z.B. des 

GutemaSes von instabilen Extrapolationsbedingungen und daher 

15 fragwurdiger Aussagekraft zeugen konnen. 

Die Abfrage der zeitlichen Konstanz einer Extrapolation kann 

bei Vorgabe eines festen GiitemaSes mit der variablen 
Extrapolationsweite (x s ) oder bei Fixierung der 

Extrapolationsweite - z. B. kann x s der Uberschreitung eines 

20 festen Grenzwertes entsprechen - liber die dann variable 
Konf idenzzahl erf olgen. 

Figur 4 zeigt einen beispielhaf ten zeitlichen Verlauf des 
GiitemaSes K(t) mit verschiedenen Werten 58. Fur diesen 
Werteverlauf wird eine Anpassungsf unktion q aus einem 
25 Kurvenvorrat ermittelt. 

Die Berechnung der jeweiligen Variablen erfolgt iterativ z.B. 
mit dem Newton-Verf ahren . 

Die zeitliche Varianziiberwachung erfolgt mit Hilfe eines 
Regressionsf its (bspw. durch ein Polynom 1 oder hoherer 
30 Ordnung, wobei die Ordnung abhangig ist von n) an die jeweils 
letzten n (z.B. n=10) zeitlichen Werte von 

Extrapolationsweiten bzw. Konf idenzzahlen. Auf diese Weise 
konnen auch sich zeitlich linear verandernde Werte der zu 
testenden Variablen prazise in ihrer Variabilitat bewertet 
3 5 werden . 

Die Bewertung der zeitlichen Varianz von Extrapolationsweite 
wird wie folgt berechnet: 
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£v •(tf(f 1 )- 9 (f))2 
c = ' 1 

Q*mitt K % 

i = Anzahl der Iterationen. 

5 

Auch hier konnen Wichtungsf aktoren Yi die aktuellsten Werte 
anders/hoher bewerten als die altesten. 
Dabei entspricht q 2 mittk y 2 mittk 

Der Wertebereich von C liegt zwischen C= 0 (keine zeitliche 
10 Streuung der Variablen, alle Variablenpunkte liegen auf der 
Regressionsgeraden) und + <*> (Variablenf olge verrauscht) . 

• Verkniipfung von V, Al norTn , S und C zum GutemaS K 

15 

Das GutemaB K verknupft die Grossen V, Al norni/ S und C bspw. 
wie folgt: 

S*C 

20 

Der Wertebereich von K liegt zwischen 0 (Extrapolation 
verwerfen) und 4- oo (Extrapolation exakt) . 
Eine Norctiierung von K kann z.B. iiber eine nichtlineare 
Beziehung erfolgen: 
25 K norm = 1 -e k . 

Bspw. ist die Anlage ein Rotor einer Gas turbine. Durch die 
Ermittelung von K des Parameters Frequenz des Rotors wird 
30 f estgestellt , dass in einer bestiiranten Zeit (xs-xn) eine 
unzulassige Grenzwertwertuberschreitung stattfinden wird. 
Daher wird die Frequenz runtergeregelt . 



35 
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Figur 5 zeigt schematise*! in einem Langsschnitt eine 
5 Gas turbine 1 . 

Als ein beispielhaf tes Bauteil ist eine Gasturbine 1 
ausgewahlt flir eine Maschine, bei der Parameter wie 
Temperatur, Schwingung, elektrische Leistung oder sonstiges 

10 Parameter gemessen werden. 

Entlang einer Welle 4 sind ein Verdichter 7, eine Brennkammer 
10 und ein Turbinenteil 13 hintereinander angeordnet. 
Das Turbinenteil 13 weist einen HeiSgaskanal 16 auf . Im 
HeiSgaskanal 16 sind Gasturbinenschauf ein 20 angeordnet. 

15 Leitschauf el- und Lauf schauf elgrenze sind wechselnd 
auf einanderf olgend vorgesehen. 

Die Gasturbinenschauf ein 20 werden beispielsweise iiber eine 
kombinierte Luft- und/oder Dampf kiihlung gekiihlt. Dazu wird 
beispielsweise dem Verdichter 7 Verdichterluf t entnommen und 
20 iiber eine Luf tzufiihrung 23 den Gasturbinenschauf ein 22 
zugefiihrt. iiber eine Dampf zuf iihrung 26 wird den 
Gasturbinenschauf ein 20 beispielsweise auch Dampf zugefiihrt. 

25 In Figur 6 ist in vereinf achter Darstellung in einem 

Langsschnitt ein Ausschnitt eines Turbinenteils 13 gezeigt. 
Das Turbinenteil 13 weist eine Welle 4 auf, die sich entlang 
einer Rotationsachse 41 erstreckt. 

Weiter weist das Turbinenteil 13 entlang der Rotationsachse 
30 41 auf einanderf olgend ein Einstrombereich 49, ein 

Beschauf elungsbereich 51 sowie ein Abstrdmbereich 53 auf. 
In dem Beschauf elungsbereich 51 sind rotierbare Lauf schauf ein 
20 und feststehende Leitschauf ein 20 angeordnet. Die 
Lauf schauf ein 20 sind dabei an der Welle 4 befestigt, wahrend 
3 5 die Leitschauf el 20 an einem die Welle 4 umgehenden 
Leitschauf eltrager 47 angeordnet sind. 
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Durch die Welle 4, den Beschauf elungsbereich 51 sowie den 
Leitschaufeltrager 47 ist ein ringf ormiger Stromungskanal fur 
ein Stromungsmedium A, z.B. HeiSdampf, gebildet. 
Der zur Zufuhr des Stromungsmediums A dienende 
5 Einstrombereich 49 wird durch ein stromauf warts des 

Leitschauf eltragers 47 angeordnetes Einstromgehause 55 in 
radialer Richtung begrenzt. 

Ein Abstromgehause 57 ist stromabwarts am Leitschaufeltrager 
47 angeordnet und begrenzt den Abstrombereich 53 in radialer 

10 Richtung, dass heiSt senkrecht zur Rotationsachse 41. 
Wahrend des Betriebs der Gasturbine 1 stromt das 
Stromungsmedium A von dem Einstrombereich 49 in den 
Beschauf elungsbereich 51, wo das Stromungsmedium unter 
Expansion Arbeit verrichtet, und verlasst danach uber den 

15 Abstrombereich 53 die Gasturbine 1. Das Stromungsmedium A 
wird anschlieSend in einen dem Abstromgehause 57 
nachgeschalteten in der Figur 5 nicht naher dargestellten 
Kondensator fur eine Dampf turbine gesammelt. 

20 Beim Durchstromen des Beschauf elungsbereichs 51 entspannt 
sich das Stromungsmedium A und verrichtet Arbeit an den 
Lauf schauf ein 20, wodurch diese in Rotation versetzt werden. 

2 5 Figur 7 zeigt in perspektivischer Ansicht eine Lauf schauf el 

20, die sich entlang einer radialen Achse 60 erstreckt. 
Die Laufschaufel weist entlang der Radialachse 60 
aufeinanderf olgend einen Bef estigungsbereich 63, eine daran 
angrenzende Schauf elplattform 66 sowie ein 
30 Schauf elblattbereich 69 auf . 

Im Bef estigungsbereich 63 ist ein Schauf elf u£ 72 gebildet, 
der zur Befestigung der Schauf el 20 an der Welle 4 einer 
Gasturbine 1 dient. Der Schauf elfuS 72 ist beispielsweise als 
Hammerkopf ausgestaltet . 

3 5 Bei herkommlichen Lauf schauf ein 20 werden in alien Bereichen 

63, 66, 69 massive metallische Werkstoffe verwendet. Die 
Laufschaufel 20 kann hierbei durch ein Gussverf ahren, durch 
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ein Schmiedeverf ahren, durch ein Frasverf ahren oder 
Kombinationen daraus gefertigt sein. Dabei weist das Bauteil 
vielfach bereits unmittelbar nach der Herstellung Fehler auf . 
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Patentansprtiche 

1. Verfahren zur Bewertung der Betriebsbedingungen einer 
Maschine (1) oder Anlage, 

5 fur die mindestens ein Parameter mehrfach gemessen wird, 

urn eine Datenbasis (6) zu bilden, 

die aus Werten (xi,yi) . . . . (x n ,y n ) des Parameters besteht, 
wobei ein GiitemaS (K) einer Extrapolation auf Grundlage 
der Datenbasis (6) berechnet wird, 
10 bei dem das GutemaS (K) eine Funktion von mindestens zwei 

Variablen der Gruppe V, Al, S, C ist, 

wobei (V) ein Verhaltnis des Wertebereichs der Datenbasis 
(6) zum Extrapolationsbereich xs ist, 
der bestimmt ist durch x s > xl, x n , 
15 wobei (Al) die x-Unscharfe der Anpassungskurve (21) in x- 

Richtung ist, 

wobei (S) die Stetigkeit als ein MaS der Veranderung der 
y-Werte der Datenbasis (6) und 

(C) die zeitliche Konstanz der Extrapolation ist. 

20 

2. Verfahren nach Anspruch 1, 

dadurch gekennzeichnet, dass 

25 die Bewertung der Betriebsbedingungen genutzt wird, 

um in Abhangigkeit des GiitemaSes (K) entsprechend den 
Parameter zu beeinf lussen. 

3. Verfahren nach Anspruch 1 oder 2, 

30 dadurch gekennzeichnet, dass 

durch die Bewertung der Betriebsbedingungen die 
Betriebssicherheit der Maschine (1) oder der Anlage erhoht 
wird, indem in Abhangigkeit des Gutemafees (K) der 
3 5 Parameter entsprechend beeinf lusst wird. 
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4 . Verf ahren nach Anspruch 1 oder 2 , 

dadurch gekennzeichnet, dass 

die Bewertung der Betriebsbedingungen fur eine Optimierung 
5 des Betriebs der Maschine (1) oder der Anlage verwendet 

wird. 



5. Verf ahren nach Anspruch 3, 
10 dadurch gekennzeichnet, dass 

fur den Parameter einen Grenzwert (18) vorgegeben wird, 
und dass ein Zeitraum ermittelt wird, 
in dem keine Uberschreitung des Grenzwertes (18) des 
15 Parameters auftritt. 



6 . Verf ahren nach Anspruch 1 , 

dadurch gekennzeichnet, dass 

20 

die Auswahl der Variablen in der Weise getroffen wird, 
dass das Giitemafi (K) von der Steigung einer 
Anpassungskurve an die Datenbasis (6) unabhangig ist. 

25 

7. Verf ahren nach Anspruch 1, 

dadurch gekennzeichnet, dass 

das GiitemaS (K) normiert ist, insbesondere durch l-e~ K . 

30 



35 8. 



Verf ahren nach Anspruch 7, 
dadurch gekennzeichnet, dass 
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das GiitemaS (K) auf einen Wertebereich von 0 bis 100 % 
normiert wird. 



5 9. Verfahren nach Anspruch 1, 

dadurch gekennzeichnet, dass 

das GutemaS (K) gegeben ist durch: 



10 K = 



10. Verfahren nach Anspruch 1 oder 9, 

dadurch gekennzeichnet, dass 

15 

das Verhaltnis (V) des Wertebereichs der Datenbasis (6) 
gegeben ist durch (X n -Xi) / (X s -Xi) . 



2 0 11. Verfahren nach Anspruch 1, 

dadurch gekennzeichnet, dass 

die Datenbasis (6) in mindestens drei Segmente (45, 48, 

52) unterteilt wird; 
25 dass fur jedes Segment (45, 48, 52) aus der Datenbasis (6) 

je ein Mittelwert g lt g 2 , g3 und je eine lineare 

Anpassungsfunktion yi, y 2 , y3 (36, 39, 42) mit den 

Steigungen Ci, c 2 und c 3 berechnet wird; 

dass ein numerisches KriimmungsmaJSes p 
30 P = gi - 2*g 2 + g 3 / 

berechnet wird, 

das die aktuelle Krurnmungsrichtung des Steigungsverlauf es 
wiedergibt ; 

dass aus einem Kurvenvorrat von Kurventypen zumindest der 

3 5 Gruppe : 
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Linear e Funktion 
Po t enz f unkt ion 
Logarithmische Funktion 
Exponentielle Funktion 



f(x)=y = a 0 + a x *x 
f(x)=lny= lnao + ai*lnx 
-» f(x)=y = a 0 + ai*ln x 
-> f (x)=ln y = lnao + a i x 



iterativ jener Kurventyp der Anpa s sung s funktion ausgewahlt 
und dem Wertebereich der gesamten aktuellen Datenbasis (6) 
angepasst wird, 

wobei der aus dem Kurvenvorrat ausgewahlte Kurventyp 
folgende Bedingungen erfullen muss: 

die Krummungsrichtung der Kurve, die aus der Extrapolation 
ermittelt wird, 

muss mit der von p ubereinstimmen 



der Quotient Q k aus Zahler (= ggf . gewichteter Mittelwert 
der Abstandsquadrate zwischen Messwerten und 
Extrapolationskurve) und Nenner (= Mittelwertquadrat des 
y-Wertebereiches der Extrapolationskurve im Bereich des 
Datenfensters) muss minimal sein: 



wobei k ein Zahler der verfugbaren Kurventypen 
(Kurvenvorrat) fur die Extrapolation ist, 
insbesondere ist y 2 m itt k = t (y m ax K + ymin K )/2]2 ' 
wobei yi(xi) der Messwert an der Stelle xi i st , 
wobei fk(xi) der Funktionswert des k-ten Kurventyps fur 
die Extrapolation an der Stelle x^ i st/ 

wobei w i ein Wichtungsfaktor fur jeden einzelnen Messwert 

oder fur alle Messwerte eines Segments ist; 

so dass die Stetigkeit (S) sich wie folgt berechnet : 



und 



&=/(*) = 



X>v,*(y t .U>/ t U-)) 2 

s mitt k I— * 1 



min (i = 1 . .min) 



S — ' 
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wobei i = 1...3 eine Numerierung fur die Segmentbereiche 
ist, 

wobei Yi : Wichtungsf aktoren l...n, 
5 wobei Oi bis O3 die Steigungen der ausgewahlten k-ten 

Kurve (36, 39, 42) fur die Extrapolation an den jeweils 
halben Segmentbreiten sind, und 
wobei Ci bis C3 die Steigungen der linearen 
Segmentsanpassungen sind. 

10 



12. Verfahren nach Anspruch 1, 

dadurch gekennzeichnet, dass 

15 die x-Unscharfe wie folgt bestimmt wird: 

Auswahl einer Extrapolationsfunktion, welche sich in 
lineare Strukturen uberfiihren lasst, also 
zumindest aus der Gruppe 

Lineare Funktion -> y = ag + a l* x 

20 Potenz funktion -> In y = lnag + a^*ln x 

Logarithmische Funktion -> y = ag + a^*ln(x) 

Exponentielle Funktion -> In y = In ag + ai*x, 

ausgew&hlt wird, 

Festlegung einer Datenbasis (6), 
25 wobei die Datenbasis (6) n korrelierte x- und y-Werte 

umf asst 

Berechnung von x und y der Datenbasis (6) sowie der 
GroSe ^x t y t 

Berechnung von c = — — ( 5^JC,y, - nxy) (i=l..n) 

** n-l 



30 



Berechnung von 5 2 = - (^*,~*) 2 (i=l..n! 

x n — 1 

Berechnung von c 2 = — — (X^"^) 2 (i=l--n 

y n-l 
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Berechnung einer Steigung b = — y 

S x 

Berechnung von a = (n-1) ( s 2 y -b 2 s 2 x ) , 

Bestimmung einer Gleichung einer Regressionsgeraden 
y = y + b(x-x) , 

mit einer Konf idenzzahl Y wird eine Grosse F (c) berechnet 

f(c) = ^(l+y), 

mit F(c) und n-2 (n = Anzahl Messwerte) Freiheitsgrade 
liefert die t-Verteilung (Student-Verteilung) einen Wert 
c, 

Bestimmung von Am 




S x <J(n-W-2) ' 

wodurch sich eine Unscharfe der Steigung m ergibt: 
b-Am<m<b + Am, 

Bestimmung der Geradegleichungen (27, 21) mit den 
Steigungen b - Am, b + Am, 

Bestimmung der Schnittpunkte (I m i n , konstant) und 
(Imax/ konstant) der Geraden mit einer Parallelen (18) 
(y=konstant) , 

die einem Grenzwert (18) entspricht, 
Bestimmung entsprechender x-Werte Ima X und Imin/ 
wobei W > I min ist, 

Berechnung von AI = I^x - I m in- 
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13. Verfahren nach Anspruch 11, 

dadurch gekennzeichnet, dass 

5 

der Wertebereich der Stetigkeit (S) im Bereich 0 und + °° 
liegt . 



10 14. Verfahren nach Anspruch 1, 

dadurch gekennzeich.net, dass 

sich das GutemaS (K) mit der Zeit verandert, und 
dass eine Bewertung der zeitlichen Varianz von 
15 Extrapolationsweite (C) wie folgt berechnet wird: 



c = ' i 

^ mitt K ^ ] y i 



20 wobei i eine Anzahl der Iterationen ist, 

wobei q 2 m itt k = [ (q ma x K + 3min K > /2 ] 2 ist 
wobei Yi e: >- n Wichtungsf aktor ist. 



25 
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Zusammenf assung 

5 Bisher wurden zur Trendanalyse visuelle Extrapolationen 

eines, bspw. kritischen, Messwertverlauf s und/oder subjektive 
Bewertungen der Auswirkungen auf technische Parameter 
vorgenommen, um die Betriebssicherheit einer Maschine 
abzuschatzen. 

10 

Die optimale Anpassung der Extrapolationsfunktion an die 
Datenbasis erfolgt erf indungsgemass nicht durch einen 
maximalen Korrelationskoef f izienten r, sondern mit Hilfe 
eines GtitemaSes (K) welches von der Steigung einer 
15 Anpassungskurve an eine Datenbasis unabhangig ist. 



(Fig. 2) 
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